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ABSTRACT 

We have made observations of 98 low-Galactic-latitude pulsars to measure pulse broadening caused 
by multipath propagation through the interstellar medium. Data were collected with the 305-m Arecibo 
telescope at four radio frequencies between 430 and 2380 MHz. We used a CLEAN-based algorithm to 
deconvolve interstellar pulse broadening from the measured pulse shapes. We employed two distinct 
pulse broadening functions (PBFs): PBFi is appropriate for a thin screen of scattering material between 
the Earth and a pulsar, while PBF2 is appropriate for scattering material uniformly distributed along 
the line of sight from the Earth to a pulsar. We found that some observations were better fit by PBFi 
and some by PBF2. Pulse broadening times (t^) are derived from fits of PBFs to the data, and are 
compared with the predictions of a smoothed model of the Galactic electron distribution. Several lines 
of sight show excess broadening, which we model as clumps of high density scattering material. A global 
analysis of all available data finds that the pulse broadening scales with frequency, v, as Td oc v~ a where 
a ~ 3.9 ± 0.2. This is somewhat shallower than the value a — 4.4 expected from a Kolmogorov medium, 
but could arise if the spectrum of turbulence has an inner cutoff at ^300-800 km. A few objects follow 
particularly shallow scaling laws (the mean scaling index (a) ~ 3.1 ± 0.1 and ~ 3.8 ± 0.2 respectively 
for the case of PBFi and PBF2), which may arise from large scale refraction or from the truncation of 
scattering screens transverse to the Earth-pulsar line of sight. 

Subject headings: ISM: structure — methods: data analysis — pulsars: general — radio continuum: 
general — scattering 



1. INTRODUCTION 

1.1. Overview 

Pulsars make excellent probes of the interstellar medium 
(ISM). Observed pulse profiles are influenced by disper- 
sion, scattering, and Faraday rotation along the line-of- 
sight (LOS) from the Earth to the pulsar. Measurements 
of pulsars in similar directions at different distances can be 
used to disentangle LOS interstellar effects and to model 
the ionized content of the ISM (Taylor & Cordes 1993; 
Bhat & Gupta 2002; Cordes & Lazio 2002a,b). 

We have undertaken multifrequency pulse profile obser- 
vations using the 305-m Arecibo telescope, concentrating 
in the swath of the Galactic plane visible from Arecibo, 
at Galactic longitudes 30° < I < 75°. The Parkes Multi- 
beam Survey (e.g. Manchester et al. 2001) has discovered 
hundreds of pulsars at low Galactic latitudes, \b\ < 5°, of 
which dozens are visible from Arecibo. Many other pulsars 



in this region are known from other survey work (e.g. Hulse 
& Taylor 1975b). Because the sensitive Multibeam Survey 
employed a higher frequency, 1400 MHz, than most other 
surveys, the "multibeam pulsars" tend to be relatively dis- 
tant and highly scattered, making them particularly useful 
for ISM studies. 

Our most fundamental measurements are the set of 
pulse shapes at different radio frequencies, from which we 
estimate the pulse broadening time scales caused by scat- 
tering, Td, for the pulsars. In addition to providing input 
data to Galactic electron density models, these measure- 
ments can be used to form an empirical relation connecting 
Td with dispersion measure, which can serve as a useful 
guide in designing large-scale pulsar surveys and in un- 
derstanding the observable population of pulsars in the 
Galaxy (e.g. Bhattacharya et al. 1992; Cordes & Lazio 
2002b). 
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The paper is organized as follows. Terminology and ba- 
sic assumptions about the ISM are summarized in § 1.2. 
Details of observations and data reduction are described in 
§ 2, and our method for deconvolving pulse broadening in 
§ 3. Our results are presented in § 4 and § 5, and in later 
sections we discuss the implications of pulse-broadening 
times for the Galactic electron density models (§ 6), as 
well as the power spectrum of electron density irregulari- 
ties (§ 7). 

1.2. Terminology and Scattering Model 

Several quantities measured by radio observations of 
a pulsar are integrals of ISM properties along the LOS 
from the Earth to the pulsar. The dispersion measure, 
DM = J ds n e (s) , is the integral of electron density, n e , 
along the LOS to the pulsar at distance D. Using a Galac- 
tic model for electron density, DM is often used to estimate 
pulsar distances. We use cataloged values of DM in the 
analysis below to estimate distances. The rotation mea- 
sure, RM = n e (s) B • ds, is the LOS integral of mag- 
netic held, B, weighted by electron density. Analysis of 
RM measurements from our data will be reported in a 
future work. 

Scattering of pulsar signals depends on fluctuations in 
the electron density, Sn e . We assume that the spectral 
density of these fluctuations follows a power-law model 
with cutoffs at "inner" and "outer" scales, k and l 0l which 
arc inversely related to the corresponding wavenumbers, 
Ki and k , by li = 2n / Ki and l Q — 2ir/n . The spectral 
density is then given by (e.g. Rickett 1977): 

0, elsewhere. 

The spectral coefficient C 2 is expressed in units of m~ 20 / 3 . 
For Kolmogorov turbulence, the spectral slope is (3 = 11/3. 

Pulse broadening is quantified by a time scale, Td, char- 
acteristic of a pulse broadening function (PBF) ht to a 
measured pulse shape. The PBF is the response of the 
ISM to a delta function. The exact form of the PBF and 
its scaling with frequency depend on the spatial distri- 
bution of scattering material along the LOS and on its 
wavenumber spectrum (Williamson 1972, 1973; Cordes & 
Rickett 1998; Lambert & Rickett 1999; Cordes & Lazio 
2001; Boldyrev & Gwinn 2003). Therefore, determina- 
tion of the PBF forms a useful means for characterizing 
the underlying scattering geometry and wavenumber spec- 
trum for scattering irregularities. The PBFs used in this 
work are described in detail in §3. 

Measured pulse scattering parameters can be related to 
the scattering measure, SM = dsC%(s), which is the 
LOS integral of C 2 . For a Kolmogorov spectrum with a 
small inner scale (e.g. Rickett 1990; Cordes & Lazio 1991; 
Armstrong, Rickett & Spangler 1995), the pulse broaden- 
ing, expressed as the mean arrival time of ray bundles (see 
Cordes & Rickett 1998), is 

(T d ) ^l.lW T SM 6/5 v- 4A D, (2) 

where v is in GHz, D is in kpc, SM is in kpc m -20 / 3 , Td is 
in ms, and W T is a geometric factor that depends on the 
LOS distribution of scattering material. 

More generally, for a power-law wavenumber spectrum, 
the broadening time scale follows a power law, 

r d otv~ a , (3) 



where (e.g. Cordes, Pidwerbetsky, & Lovelace 1986; Ro- 
mani, Narayan, & Blandford 1986), 



(1) 



2/3 
(0-2) 



f3 < 4 
P > 4 



(4) 



Thus, determination of a yields information about the 
wavenumber spectrum. For a Kolmogorov spectrum, (3 = 
11/3, implying a = 4.4, and Eq. 3 reduces to 2. This re- 
sult holds if the inner scale of the spectrum is too small 
to influence the measurements (Cordes & Lazio 2002). As 
we discuss later, we infer that the inner scale likely does 
influence some of the scattering measurements. 

Finally, the decorrelation bandwidth, 1/4, is related to Td 
by 2-KT d Vd = Ci, where the constant C\, of order unity, 
depends on the geometry of the scattering material and 
the form of the wavenumber spectrum (Cordes & Rickett 
1998). 

2. OBSERVATIONS AND DATA REDUCTION 

The observations were made at the Arecibo Observa- 
tory. New data for 81 pulsars were obtained in several 
observing sessions from 2001 May to 2002 November. For 
the analysis in this paper we also use the data collected 
by Lorimer et al. (2002) for 17 pulsars, yielding a total of 
98 pulsars. We concentrated on pulsars for which pulse 
broadening observations had not previously been made. 
Prominent among these are 38 discovered in the Parkes 
Multibeam Survey (e.g. Manchester et al. 2001), 30 from 
the Hulse- Taylor survey, including 17 with new timing 
solutions (Lorimer et al. 2002), and 30 others (Taylor, 
Manchester, & Lyne 1993; Hobbs & Manchester 2003) 

Data acquisition systems used for the observations are 
summarized in Table 1. Signals were collected separately 
at four radio frequencies, 430, 1175, 1475 and 2380 MHz. 
The range of frequencies was chosen to allow detection 
of pulse broadening over a wide variety of pulsar scatter- 
ing measures; specific frequencies were chosen according to 
receiver availability and radio frequency interference envi- 
ronment. The strong dependence of Td on frequency im- 
plies that, for most objects, pulse broadening will be mea- 
surable at only a subset of the four frequencies. For pulsars 
with little scattering, pulse broadening is detectable only 
at the lowest frequency, if at all. By contrast, for pulsars 
with heavy scattering, broadening may be measurable at 
high frequencies and may be so large as to render pulsa- 
tions undetectable at lower frequencies. When pulsations 
are undetectable at 430 MHz, the cause may also involve 
a combination of relatively small flux density and large 
background temperature. 

In the absence of any prior knowledge of flux density at 
the higher frequencies, we adopted fixed integration times 
for all objects in a first pass of observations. Based on 
the initial results, one or more re-observations were made 
during later sessions for those objects and frequencies with 
low signal-to-noise ratios. 

Observations at 430 MHz were made with the Penn 
State Pulsar Machine (PSPM), an analog interbank spec- 
trometer providing 128 spectral channels spanning an 8 
MHz band in each of two circularly polarized signals. 
Power measurements in each channel were synchronously 
averaged in real time at the topocentric pulse period, yield- 
ing pulse profiles with time resolution of approximately 1 
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milliperiod. Dedispersion was done off line, reducing each 
observation to a single pulse profile. 

Observations at 1175, 1475, and 2380 MHz were made 
with the Wideband Arecibo Pulsar Processor (WAPP), a 
fast-dump digital correlator (Dowd, Sisk, & Hagen 2000). 
Input signals to the WAPP were digitized into three lev- 
els and output correlations were accumulated and written 
to disk as 16-bit integers. We recorded long time series 
of auto-correlation functions (ACFs) and cross-correlation 
functions (CCFs) of the two circularly-polarized polariza- 
tion channels. The ACFs were used in the analysis of this 
paper. A polarization analysis that utilizes both the ACFs 
and CCFs will be reported in a subsequent paper. 

In off-line analysis, the ACFs were van Vleck corrected 
(e.g. Hagen & Farley 1973), Fourier transformed, dedis- 
persed, and synchronously averaged to form average pulse 
profiles with, typically, 1 milliperiod resolution. The soft- 
ware tools used to analyze the WAPP and PSPM data are 
described by Lorimer (2001). 

Figure 1 shows the pulse profiles obtained from our mul- 
tiple frequency data. Profiles for 5 pulsars are not shown 
due to poor data quality. For the 37 multibeam pulsars 
shown, these represent the first observations at frequencies 
other than 1400 MHz, and, in almost all cases, signal-to- 
noise ratios for the profiles are superior to those obtained 
in the original multibeam survey data. For nearly all 39 
previously known pulsars in the sample shown here, the 
profiles are the best quality profiles obtained to date. 

3. DECONVOLUTION METHOD 

We used a CLEAN-bascd method (Bhat et al. 2003) 
for deconvolving scattering-induced pulse-broadening from 
the measured pulse shapes. This method does not rely on 
a priori knowledge of the pulse shape, and it can recover 
details of the pulse shape on time scales smaller than the 
width of the PBF. A number of trial PBFs may be used, 
with varying shapes and broadening times, correspond- 
ing to different LOS distributions of scattering material. 
The "best fit" PBF and broadening time are determined 
by a set of figures of merit, defined in terms of positivity 
and symmetry of the final CLEANed pulse, along with the 
mean and rms of the residual off-pulse regions. Details of 
the method and tests of its accuracy are given in Bhat et 
al. (2003). 

We used two trial PBFs. The first, PBFi, is appropriate 
for a thin slab scattering screen of infinite transverse ex- 
tent within which density irregularities follow a square-law 
structure function 1 (Lambert & Rickett 2000). The PBF is 
given by a one-sided exponential (Williamson 1972, 1973), 

PBFi(i) = t^ 1 exp(-t/T d )U(t), (5) 

where U(t) is the unit step function, U(< 0) = 0, U(> 
0) = 1. This function has been commonly used in previous 
pulsar scattering work. 

The second broadening function, PBF2, corresponds to 
a uniformly distributed medium with a square-law struc- 
ture function. This PBF has a finite rise time and slower 
decay, 

PBF 2 (t) = (7r 5 r//4t 5 ) 1/2 cxp(-7r 2 r d /At)U{t). (6) 

1 

The spatial structure function Dp(s) of a quantity F(x) is defined as 
Dp(s) = {(F(x + s) — F(x)) 2 ) where s is the spatial separation (lag 
value) . 



This PBF is a generic proxy for more realistic distributions 
of scattering material. 

Additional PBFs, not used in our analysis, include those 
for media with Kolmogorov wavenumber spectra, which 
can yield non-square-law structure functions (e.g. Lambert 
& Rickett 1999), and scattering screens that are truncated 
in directions transverse to the LOS, as may be the case for 
filamentary or sheet-like structures, which have PBFs that 
correspondingly are truncated at large time scales (Cordes 
& Lazio 2001). 

Note that the pulse broadening time, Td, has different 
meanings for PBFi and PBF 2 . For PBFi, Td is both the 
e _1 point of the distribution and the expectation value of t. 
For PBF 2 , Td is close to the maximum of the distribution, 
which is at (n 2 /10)Td = 0.997^, while the expectation value 
of t is (7r 2 /2)r d = 4.93r d . 

For some of the pulsars we obtained an acceptable fit 
to Td using both PBFi and PBF 2 , while in others only 
one of the PBFs provided an acceptable fit. Acceptable 
fits were those that yielded deconvolved pulse shapes that 
were positive, semi-definite; we reject cases which yielded 
unphysical pulse shapes (such as profiles with negative go- 
ing components). In many cases the pulse broadening is 
not large enough to be measured, in which case we quote 
upper limits on T d (see Table 2). 

As noted earlier in this section, our method relies on a 
set of figures of merit for the determination of the best fit 
PBF for a given choice of the PBF form (see Bhat et al. 
(2003) for details). Among the different parameters used 
for this determination, the parameter f r is a measure of 
positivity, and can serve as a useful indicator of "good- 
ness" of the CLEAN subtraction. However we emphasize 
that the absolute value of this parameter may also depend 
on the degree of scattering, the noise in the data, shape 
of the intrinsic pulse, etc., and therefore a comparison of 
the results for different data-sets will not be meaningful. 
Nonetheless it can still be used for a relative comparison of 
the results obtained using different PBFs for a given pulse 
profile. For successfully deconvolved pulses, we expect 
f r ~l; larger values imply slightly ovcrCLEANed pulses. 
Based on this approach, the PBF with a lower value of f r 
can be considered to be the better of the two PBFs. 

4. DERIVED INTRINSIC PULSE SHAPES 

Figures 2 and 3 show results from the CLEAN-based de- 
convolution of our data. In each panel, the best-fit PBF is 
shown along with the observed pulse shape and the decon- 
volved (intrinsic) pulse shape. As is evident in the figures, 
the derived pulse shapes are much narrower and signif- 
icantly larger in amplitude than the observed ones. In 
several cases, the deconvolved pulse shapes reveal signifi- 
cant structure which is not easily visible in the measured 
pulse profiles. For example, PSRs J1913+0832 (Fig. 2) 
and J1858+0215 (Fig. 3) at 1175 MHz, have derived pulse 
shapes that are distinct doubles, a property that is al- 
most entirely masked by broadening in the raw profiles. 
In several other cases (e.g. PSR J1912+0828 at 1175 
MHz, Fig. 2; PSR J1927+1852 at 430 MHz, Fig. 3), the 
measured pulse shapes show faint signatures of a double, 
which are confirmed and reinforced by the deconvolution 
process. Data for PSR J1906+0641 at 1175 MHz (Fig. 2) 
and PSR J1942+1743 at 430 MHz (Fig. 3) show that the 
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technique yields details of complex, multi-component pulse 
shapes. 

While the deconvolution algorithm usually produces ac- 
curate pulsar profiles, some cautions are in order. As we 
already discussed, for many objects successful deconvolu- 
tion is possible using both PBF forms. Figures 2 and 3 
include several examples of this kind. In some cases, sig- 
nificantly different intrinsic pulse shapes result from de- 
convolution with the two different PBFs. The data for 
PSR J1852+0031 show an extreme example of this: de- 
convolution with PBFi yields a pulse shape with three 
merged components (Fig. 2), while use of PBF 2 yields a 
distinct double pulse shape (Fig. 3). Further examples of 
substantial discrepancies include PSR J1858+0215 at 1175 
MHz, where use of PBF2 yields a double pulse, while use 
of PBFi yields a simpler, single-peaked pulse profile, and 
PSR J1853+0545 at 1175 MHz, where use of PBF 2 yields 
a much narrower and more featureless profile than PBFi. 
However, such cases are exceptions rather than the rule. 
There are many examples where nearly identical intrin- 
sic pulse shapes result with either of the two PBFs. Data 
from PSRs J1905+0616 and J1907+0740 at 430 MHz, and 
J1908+0839 and J1916+1030 at 1175 MHz, belong to this 
category. 

We emphasize that we have used only two extreme ex- 
amples from the infinite set of possible PBFs, and there 
may be LOSs for which other PBFs would be more appro- 
priate. An exhaustive analysis using additional PBFs is 
beyond the scope of this paper, though such an analysis 
may be valuable in identifying important aspects of the 
ionized interstellar medium. In the remainder of the pa- 
per, we focus on measurements of pulse-broadening times 
and their implications for models of Galactic free electron 
density. 

5. PULSE-BROADENING TIMES 

Our estimates of Td are summarized in Table 2. The 
columns are: (1) pulsar name, (2) reference, (3) pulse pe- 
riod, (4) DM, (5) Galactic longitude, (6) Galactic latitude, 
(7) observation frequency, (8) estimate of Td using PBFi, 
and (9) its figure of merit (/ r ), (10) estimate of T<j using 
PBF 2 , and (11) its figure of merit (/ r ), (12) model estimate 
of pulse broadening using the TC93 model (^^93), and 
(13) model estimate of pulse broadening using the NE2001 
model (T(2 jne 2ooi)- The definitions of model estimates are 
discussed below along with a comparison with measured 
values of Td- We successfully measured Td for 56 of the 98 
target objects, (of which 15 have measurements at more 
than one frequency), and obtained upper limits on Td for 
31 objects. 

5.1. Scaling of Td with Frequency 

For 15 pulsars, we have measured Td at more than one 
frequency, typically at 1175 and 1475 MHz, but in one 
case, PSR J1853+0545, also at 2380 MHz. For 12 of these, 
estimates of Td were possible using both PBFi and PBF 2 
(Table 2), and we derive the estimates of the scaling index 
a in both cases (columns 4 and 6 in Table 3). We use 
ai and a 2 to denote the scaling indices for the two PBF 
cases, PBFi and PBF 2 , respectively, and the correspond- 
ing values for (3 (obtained by use of Eq. 4) are denoted as 
Pi and /3 2 . For PSR J1853+0545, measurements of Td are 



available for 3 frequencies, and we estimate a for all three 
pairs of frequencies. 

Despite the small sample of measurements, we find: (1) 
most cases show significant departures from the traditional 
v~ iA scaling expected for Td, and (2) the inferred scaling 
index depends on the type of the PBF adopted for the 
deconvolution. These have important implications for the 
interpretations that will ensue in terms of the nature of 
the wavenumber spectrum, as we discuss below. 

5.2. Scaling of Td with DM 

An empirical relation connecting the pulse-broadening 
time and DM serves as a useful guide in designing large- 
scale pulsar surveys. An ideal pulsar survey will be scatter- 
ing limited rather than dispersion limited. Most surveys 
to date have not, in fact, been scattering limited; this is 
not because they have been poorly designed, but rather 
because they have been constrained by data throughput 
and computational limitations. In other words, scatter- 
ing plays a significant role in determining the maximum 
distance to which a pulsar can be detected and thus influ- 
ences the observable population of pulsars. The relation 
also provides some useful insights into the large-scale dis- 
tribution of free electrons (n e ) and the strength of their 
density fluctuations (5n e ) in the Galaxy. 

Figure 4 shows a scatter plot of Td and DM. Most of the 
points at smaller DMs (s 100 pc cm~ 3 ) are derived from 
measurements of dccorrclation bandwidth, i/<j, which are 
converted to scattering times by Td — C\/2-KVd, assuming 
C\ = 1. Direct measurements of pulse- broadening times 
dominate at larger DMs (^ 100 pc cm -3 ). Evidently, there 
is a strong correlation between DM and Td over the 10 or- 
ders of magnitude of variation in r^. The values of DM 
cover only 3 orders of magnitude, signifying a strong de- 
pendence of pulse-broadening on DM. There is also large 
scatter of Td about the trend, roughly 2 to 3 orders of 
magnitude. Some of the scatter results from the fact that 
we have scaled all measurements to a common frequency 
of 1 GHz using Td oc v~ iA . However, alternative scal- 
ing indices will yield an error of no more than about 0.4 
in log Td- At lower DMs, some of this scatter may be at- 
tributed to refractive scintillation effects which cause fluc- 
tuations in the decorrelation bandwidth (e.g. Bhat, Rao, 
& Gupta 1999a). Also, some of the scatter may be due 
to the small numbers of "scintles" contributing to a mea- 
surement. At larger DMs, the scatter is primarily caused 
by strong spatial variations in C 2 . 

We fit the values of Td and DM shown in Figure 4 using 
a simple parabolic curve of the form (e.g. Cordes & Lazio 
2002b) 

log T d ~a + 6(log DM) + c(log DM) 2 - a log u, (7) 

where v is the frequency of observation in GHz, and Td 
is in ms. Previous work has assumed a fixed scaling in- 
dex, a = 4.4, while fitting for the coefficients a, b and 
c. In the light of our results discussed in § 5.1, and also 
other recent work (e.g. Lohmer et al. 2001) that suggest 
a departure from the traditional r<j oc v~ 4A behavior, we 
treat the scaling index a as an additional parameter in 
determining the best fit curve. Note that most published 
measurements of Td were determined by assuming a PBF of 
the form PBFi (and assuming the conventional frequency- 
extrapolation approach). Hence we use our values of Td 
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determined by using the same form of PBF (column 8 of 
Table 2) in order to ensure uniformity of the data used for 
the fit. Furthermore, to allow an unbiased fit for a, we use 
measurements in their unsealed form, i.e., direct estimates 
of Td and Vd at the observing frequencies 2 . The data used 
for our fit, many from prior compilations, include 148 es- 
timates of Vd and 223 estimates of Td (of which 64 are our 
own measurements), thus 371 measurements in total. Note 
that the upper limits are excluded from the fit, as none of 
them seem to impose any constraints to the fit. For a sub- 
set of these objects, measurements are available at multiple 
frequencies. The best-fit curve from our analysis is shown 
as the solid line in Figure 4. Our re-derived coefficients, 
a = -6.46, b = 0.154 and c = 1.07, arc only slightly dif- 
ferent from those of Cordes & Lazio (2002b), a — —6.59, 
b = 0.129 and c = 1.02. Interestingly, the global scaling 
index derived from our best fit is a — 3.86 ± 0.16, which is 
significantly less than the canonical value of 4.4 appropri- 
ate for a Kolmogorov medium with negligible inner scale. 

There are several plausible explanations for departure 
from the v~ 4A scaling behavior for Td, such as (i) the 
presence of a finite wavenumber cutoff associated with an 
inner scale, (ii) a non-Kolmogorov form for the density 
spectrum, and (hi) truncation of the scattering medium 
transverse to the LOS, as addressed by Cordes & Lazio 
(2001). Presently available observational data suggest that 
option (i) may apply, so we investigate the effects of an 
inner scale on the scaling laws for r^. Option (iii) may 
be relevant for specific LOSs that contain filamentary or 
sheet-like structures that could mimic truncated screens. 
In addition, there is yet another effect whereby a weaken- 
ing of the scaling index (as deduced from measurements of 
Vd and Td) could result from refraction effects in the ISM. 
As argued theoretically and demonstrated through obser- 
vational data, refraction from scales larger than those re- 
sponsible for diffraction will bias the diffraction bandwidth 
downward, corresponding to an upward bias on pulse- 
broadening (e.g. Cordes, Pidwerbetsky, & Lovelace 1986; 
Gupta, Rickett, & Lyne 1994; Bhat et al. 1999b). The 
refraction effects will be stronger at higher frequencies as 
one approaches the transition regime between weak and 
strong scattering, which will tend to weaken the frequency 
dependence from 4.4 to a lower index. For pulsars at low 
DMs (say, DM S 100 pc cm~ 3 ), this transition is expected 
near ~l-3 GHz. Our sample contains many low-DM ob- 
jects with measurements at ~l-2 GHz where such an effect 
may be significant. 

5.2.1. Effect of finite inner scale on a 

The presence of a finite inner scale can potentially mod- 
ify the frequency scaling index as estimated from measure- 
ments of Vd and Td- Cordes & Lazio (2002b) show that 
these effects become apparent above a "crossover point" 
that is a function of distance (or DM) as well as the ob- 
serving frequency v. The crossover point can be defined 
for commonly used observables such as 9d (angular broad- 
ening), Vd and Td- In order to examine our data for any 
such signatures of an inner scale, we define a "test quan- 
tity" in terms of Td, v, and distance (D) that is directly 
related to the inner scale, expressed in units of 100 km, 

2 

Note that many published data, such as those in Taylor et al. (1995), 
are already pre-scaled to a common frequency of 1 GHz. 



'ioo = £»/(100km). The crossover point Td iCross is related 
to the inner scale by (see Eq. A20 of Cordes & Lazio 
(2002b)), 

Td, cross « 5.46 ms -D(^ghz^ioo)~ 2 , (8) 
where D and v are in kpc and GHz, respectively. Thus, 
a useful test quantity for identifying a break point in the 
frequency scaling is Td, cross v 2 /D. In the analysis that 
follows, we use a simple linear relation to convert DM 
measurements to distances, D — DM/(1000(ra e )), where 
(n e )=0.03 cm -3 is the mean electron density, and DM is 
in units of pc cm~ 3 . We emphasize that we adopt such a 
simplistic approach as a preliminary step, and will defer 
to another paper a more detailed and complete analysis 
using proper electron density models and the independent 
pulsar distance estimates. 

We split the data-set into two parts, below and above a 
chosen break point value for this test quantity, and for each 
case we refit the parabolic curve in Eq. 7 for the best fit a 
while keeping the coefficients a, 6, c fixed at their global fit 
values. We do this exercise for several break point values 
in the range 0.03-3.3, determining the difference in best fit 
a values for the two samples in each case {5a = aw — abh , 
where uu and abh denote the values of a for the samples 
that are below and above the break point). If an inner 
scale effect is truly relevant, we will expect a significant 
difference in a for the two samples (with a larger value for 
the sample below the break point, i.e., au > abh )■ 

Figure 5 shows a plot of da vs. the test quantity 
Td, cross ^ 2 / D , along with a corresponding plot of the best fit 
X 2 (x 2 = Xi +xl ) where Xi and \\ denote the correspond- 
ing values of \ 2 f° r t ne two data-sets). The maximum in 
8a roughly coincides with the minimum in x 2 , suggest- 
ing that the inner scale effect is real. Our analysis shows 
a sharp minimum for \ 2 at log(rd. cross ^ 2 /D) « —0.57 
(Fig. 5). Formally, the ±lcr error on the break point value 
of log( 

T d,cross^ 2 /-D) is ±0.05. However, the valley in \ 2 is 
much broader than implied by this error. We take a more 
realistic range to be —1 to —0.3 in the log, corresponding 

to an inner scale k « 100 km (5.46L>/r d ^ 2 ) 1/2 w 300 to 
800 km. 

The broadness of \ 2 is caused in part by our assump- 
tion of a simple proportionality between distance and DM 
and also by the likely variation of inner scale between lo- 
cations in the Galaxy. Some theories for density fluctu- 
ations in the ISM would associate the inner scale with 
the proton gyroradius for thermal gas. The gyroradius is 
r g « 1658 kmT 4 1/2 B"^ for a temperature T = 10 4 T 4 K and 
a magnetic field strength B expressed in micro gauss. For 
ionized gas in the warm phase of the ISM, we expect the 
temperature to vary by a factor of 2 to 4 and the field 
strength by at least a similar factor. Thus, we would ex- 
pect the gyroradius to vary by at least a factor of five, 
which is not inconsistent with the appearance of \ 2 m 
Figure 5. Given the expected variation of the gyroradius 
in the ISM, it is perhaps surprising that we see any kind 
of minimum in \ 2 at all. 

Several authors have investigated the effect of an inner 
scale, and constraints are available from various kinds of 
observations. For example, Spangler & Gwinn (1990) de- 
rived an inner scale of ~50-200 km from an analysis of in- 
terferometer visibility measurements from VLBI observa- 
tions, which are, interestingly, of the order of our estimates 
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derived from pulse-broadening data. Further, as noted by 
Moran et al. (1990), observations of NGC 6334B (the ob- 
ject with the largest known scattering disk) at centimeter 
wavelengths are consistent with an inner scale larger than 
35 km. Studies of long-term flux density variations of pul- 
sars at low radio frequencies, however, indicate a much 
larger inner scale (e.g. ~ 10 2 -10 4 km from the work of 
Gupta, Rickett, & Coles (1993)). While some discrepan- 
cies prevail between the estimates deduced from different 
observations, it appears that effects due to an inner scale 
are well supported by a number of observations. 

6. GALACTIC ELECTRON DENSITY MODELS 

Our sample largely comprises high-DM, distant pulsars 
and hence provides useful data for improving upon electron 
density models for the inner parts of the Galaxy. We com- 
pare our data with predictions from both the TC93 (Tay- 
lor & Cordes 1993) and NE2001 (Cordes & Lazio 2002a,b) 
models, which yield values for SM that may be used in 
Eq. 2. The newer model, NE2001, has made use of only 
DM values of some of the multibeam pulsars; hence, our 
measurements of Td allow an independent test of the new 
model. 

Figures 6 and 7 show plots of the measurements of pulse- 
broadening times against the predictions from the new 
and old electron density models, respectively. In order 
to examine more general trends, we also plot all the pub- 
lished measurements (see Taylor et al. (1995) and refer- 
ences therein), after scaling to a common frequency of 1 
GHz using Td oc v~ iA . A significant number of measure- 
ments show reasonable agreement with the model predic- 
tions, suggesting that the models depict fairly good rep- 
resentation of the large-scale picture in the Galaxy. How- 
ever, significant discrepancies are evident in many cases, 
compared against the predictions from either of the two 
models. For a majority of the measurements (~75%) from 
our own observations, the discrepancy is significantly lower 
with the predictions of NE2001 than with those of TC93. 
In some cases, the agreement with the model prediction 
shows improvements of the order a factor of two or better. 
Given that our measurements were not part of the inputs 
for the new model, this comparison makes an independent 
test of the new model. 

6.1. Clumps of Excess Scattering 

As discussed in § 1.2, the measured broadening time Td 
is related to the total amount of scattering, usually quanti- 
fied as the scattering measure, SM (see Eq. 2). For a given 
scattering geometry (indicated by the corresponding geo- 
metric factor W T in the equation) , we can invert this equa- 
tion to derive the scattering measure. We assume W T =1 
and estimate the effective SM for a uniform medium. The 
estimated values of SM (SM moas , in the conventional units 
of kpc m~ 20 / 3 ) are listed in Table 4 (column 7). 

Figure 8 shows the distribution of inferred SMs at the 
locations of pulsars. The spiral arm locations are adopted 
from the NE2001 model and the pulsar distances are the 
revised estimates using this new electron density model. 
A more useful quantity is the departure of the mea- 
sured quantity (jd or SM) from the model predictions. 
In order to examine this in detail, we plot the quantity 
| log(r t ;/r ( j in e200i)| at the locations of the pulsars (Fig. 9). 



Significant departures are seen toward many LOSs. In 
the case of low-DM pulsars, these may be due largely to 
measurement errors due to refractive scintillation effects 
(Bhat, Rao, & Gupta 1999a). For distant, high-DM pul- 
sars, departures from the model predictions are in general 
larger in the inter-arm region. Most published data are in 
good agreement with the model predictions as expected, 
while several of the new measurements differ significantly 
from the model expectations. 

A closer examination of Figures 6 and 7 reveals that 
despite the general agreement seen with a large number 
of measurements, the models still underestimate the to- 
tal amount of scattering for many LOSs. The underes- 
timates are accounted for easily by relaxing one or more 
assumptions that underly the calculation of SM and its 
interpretation, as has been pointed out by Cordes & Lazio 
(2002b). In particular, clumps of enhanced scattering are 
likely due to unmodeled features associated with Hn re- 
gions or supernova shocks. Following Cordes & Lazio 
(2002b) (see also Chatterjee et al. (2001)), we characterize 
such "clumps" in terms of the incremental SM and DM 
due to them. We calculate the increments associated with 
a clump as 



<5DM = n ec Ss 



and 



SSM = C 2 c Ss 



(9) 



(10) 



where n ejC is the mean electron density and C 2 c is a mea- 
sure of the fluctuating electron density inside a clump and 
Ss is the size of the clump region. The parameter C 2 c can 
be expressed in terms of the electron density and the 
"fluctuation parameter," F c (see TC93; Cordes & Lazio 
(2002a)), 

Cl c =C SM F c nl c (11) 

where Csm is a numerical constant that depends on the 
slope of the wavenumber spectrum, and is defined as 

Csm — [3(27r) 1 / 3 ] K u for a Kolmogorov spectrum, 
where the scale factor K u = 10.2 m~ 20 / 3 cm 6 yields SM 
in the conventional units of kpc m~ 20 / 3 . The fluctuation 
parameter F c depends on the outer scale (Z Q ), filling fac- 
tor (77), and the fractional rms electron density inside the 
clump. It is defined as (TC93) 

F c = CeVX" 273 , (12) 

where £ — (n e 2 )/(n^) 2 , and e = ((Sn e ) 2 ) /n^ 2 . From equa- 
tions 9-11, the ratio of increments in SM and DM is given 

by 

^=Cs M F c n eiC . (13) 

The above expression can be re-written as 

F c (<5DM) 2 



SSM = C SM - 



Ss 



(14) 



For large distances of a few to several kpc that are relevant 
for our measurements, it is a fair assumption that the LOS 
to the pulsar may encounter several such clumps. Assum- 
ing a clump thickness of ~10 pc (typical size of known H11 
regions), and a volume number density for clumps, n c ~1 
kpc -3 , we obtain the values of F C (5DM) 2 for the subset 
of measurements in Table 4 that show excess SM (see also 
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Fig. 10). The constraints derived from our data lie within 
a broad range of 

10 -5.3 < f^sbM) 2 < lO- 18 , (15) 

which is consistent with values needed to account for the 
excess scattering toward the LOS to pulsar B09 19+06 de- 
rived by Chatterjee et al. (2001). If we assume a fluctu- 
ation parameter (F c ) of 10 for the clumps, which is con- 
sistent with values in TC93 and Cordes & Lazio (2002a), 
the required range in 5DM is 7 x 10~ 4 < <5DM < 4 x 10~ 2 
pc cm~ 3 . For the assumed size of 10 pc for the clumps, 
this implies 10~ 5 ;S n e c s 4x 10~ 3 cm~ 3 . In reality, both the 
fluctuation parameter F c , as well as the sizes (5 s) and num- 
ber of clumps (n c ), will vary with the LOS; nonetheless, 
the inferred values of <5SM arc such that the derived con- 
straints on the clumps are well within the range of physical 
possibilities. Note also that the implied perturbations of 
DM are rather small, a fact that highlights the situation 
that relatively small changes in the local mean electron 
density can translate into large changes in the amount of 
scattering. 

7. IMPLICATIONS FOR THE ELECTRON DENSITY 
WAVENUMBER SPECTRUM 

Our measurements of the scaling index a and the im- 
plied power spectral slopes (3 are summarized in Ta- 
ble 3. In a few cases, the estimates of a are consistent 
with the simple, Kolmogorov scaling of a = 4.4 (e.g. 
PSRs J1853+0545, J1913+1145, and J1920+1110). How- 
ever, in most cases the measured scaling is significantly 
weaker than even I/" 4 (e.g. PSR J1856+0404). Overall 
the measurements show a possible departure from the tra- 
ditional expectation of v~ scaling, with a mean scaling 
index (a) « 3.12 ± 0.13 using the results for PBFi, and 
« 3.83 ± 0.19 using those for PBF 2 , in agreement with 
other recent work (Lohmer et al. 2001), and also compa- 
rable to a global scaling index 3.86 ± 0.16 inferred from 
our parabolic fit to Td vs DM data. 

Figure 11 summarizes the current state of the estimates 
of a derived from measurements of decorrelation band- 
widths and pulse-broadening times. In addition to the 
present work which yielded a for 15 objects, this includes 
the recent measurements from Lohmer et al. (2001) (for 
9 high-DM objects) and those from Cordes, Weisberg, & 
Boriakoff (1985) (for 5 objects at low DMs) derived from 
measurements of Ud- Barring a few outlier cases, it appears 
that the scaling index is lower for objects of larger DMs 
(£ 200 pc cm -3 ), while it seems consistent with the Kol- 
mogorov expectation for objects at lower DMs (although 
these are only 5 in number). A similar result is also in- 
dicated by our analysis of the DM dependence of pulse 
broadening times (§ 5.2.1). 

We now return our attention to the dependence of the 
scaling index on the PBF form adopted for the analysis. 
PSR J1853+0545 is a particularly illustrative example. 
For this object, based on the results for PBFi, we estimate 
a mean scaling index much lower than 4.4, (a) = 3.1 ±0.2. 
However, use of PBF 2 yields scaling indices that agree well 
with that expected for a j3 = 11/3 spectrum. Naturally, 
the two cases may lead to widely different interpretations 
in terms of the nature of the wavenumber spectrum. Sim- 
ilarly, for PSR J1857+0526, while the estimate of a de- 
duced from Td values obtained for the case of PBFi imply 



a power-law index that approaches the Kolmogorov value, 
even if it is a little low, the results for the case of PBF 2 
yield a much lower value. 

Given all of this, it is important to attempt to use an 
approximately correct form for the PBF before attempt- 
ing any serious interpretation in terms of the nature of the 
spectrum. A mere departure from the expected v scal- 
ing need not necessarily signify an anomaly for the scat- 
tering along that LOS. However, in general it is clearly 
difficult to know what is the correct form of the PBF for a 
given LOS. Additional figures of merit such as the derived 
intrinsic pulse shapes, in particular their dependence with 
frequency, and the number of CLEAN components (Bhat 
et al. 2003) may help to resolve ambiguities in some cases, 
but the general problem remains a difficult one. 



8. SUMMARY AND CONCLUSIONS 

We have used multi-frequency radio data obtained with 
the Arecibo telescope for a sample of 98, mostly dis- 
tant, high-DM, pulsars to measure in particular the pulse- 
broadening effect due to propagation in the inhomoge- 
neous interstellar medium. For 81 of these objects we 
obtained data at 0.4, 1.2, 1.5 and 2.4 GHz, while data 
for the remaining 17, at 0.4, 1.2 and 1.5 GHz, are from 
the recent work of Lorimer et al. (2002). We employed 
a CLEAN-based deconvolution method to measure pulse- 
broadening times. In this process we tested two possible 
forms of the pulse-broadening function that characterizes 
scattering along the LOS. As a by-product, the method 
also yields estimated shapes of the intrinsic pulse profiles. 

The present work has resulted in new measurements of 
pulse-broadening time for 56 pulsars, and upper limits for 
31 pulsars. These data, along with similar measurements 
from other published work, were compared with the pre- 
dictions from models for the Galactic free electron density. 
New measurements allow an independent test of the elec- 
tron density model recently developed by Cordes & Lazio 
(2002a). While a majority of the data is in reasonable 
agreement with the model predictions, evidence for excess 
scattering is seen for many LOSs. We consider the possi- 
bility whereby the excess scattering can be accounted for 
by using "clumps," regions of enhanced scattering in the 
Galaxy. Depending on the distance, a given LOS may 
contain one or more of such clumps, and we derive useful 
constraints on their properties. 

For a small subset of objects, our data also allow es- 
timation of the frequency scaling indices for the pulse- 
broadening times, most of which show significant depar- 
tures from the traditional v~ AA behavior expected for the 
case of a Kolmogorov power-law form for the spectrum of 
density irregularities. Our analysis also suggests that the 
inferred scaling indices depend on the type of PBF adopted 
for the analysis. We combined our data with those from 
published work to revise the empirical relation connecting 
pulse-broadening times and dispersion measures. The in- 
ferred frequency scaling index from a globat fit is 3.9±0.2, 
less than that expected for the case of a Kolmogorov spec- 
trum. Our analysis also suggests the possibility of an inner 
scale in the range ^300-800 km for the spectrum of tur- 
bulence. Further, the intrinsic pulse shapes deduced from 
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our analysis for several of the pulsars are likely to be com- 
parable to the actual pulse shapes, and hence may prove 
useful for applications such as the study of pulsar emission 
properties. 
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APPENDIX 



APPENDIX A 
PROFILE DATABASE 

The basic data (i.e. the pulse profiles in Fig. 1) presented in this paper are also available as an electronic data set. 
The full database is packaged as a gzipped tar file, AOmultif req_prof s . tar . gz (which includes 345 pulse profiles from 
our observations), and is available over the Internet http://web.haystack.mit.edu/staff/rbhat/aoprofs or can be 
downloaded via anonymous ftp from the ftp site web.haystack.mit.edu (the directory is pub/rbhat/aoprof s). These 
profiles are stored as individual files in simple ascii format, which consists of a header line of the basic observing parameters 
followed by an ascii list of pulse bin number and the intensity value (in arbitrary units) in a two-column format. Each 
file is given a generic name of the format pulsar .freq.prf, where pulsar is the name of the pulsar and freq is the 
frequency of observation in MHz. An example header is shown below, along with a description of the various parameters 
included in the header. 

# mjdobs mjdsec per np freq refdm nbin siteid scanid source 



where 

mjdobs : Date of observation (MJD) 

mjdsec : Time of observation (seconds, UTC) with 

respect to mjdobs 

per : Pulse period (seconds) 

np : Pulse count 

freq : Frequency of observation (MHz) 

refdm : Dispersion measure (pc cm~ 3 ) 

nbin : Number of bins in the pulse profile 

siteid : Site ID of observations ('3' for Arecibo) 

scanid : Scan number of observation 

source : Source name 
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Table 1 
Data acquisition parameters. 



Frequency 


Bandwidth 


Time 


Spectral 


Instrument 


Integration Timc a 


(MHz) ' 


(MHz) 


Resolution 


Channels 




per Scan (min) 


430 


8 


10~ 3 P 


128 


PSPM 


10 


1175 


100 


256 /is 


256 


WAPP 


5 


1475 


100 


256 ^s 


256 


WAPP 


5 


2380 


50 


256 /us 


64 


WAPP 


10 



a Multiple scans were made for pulsars with low signal-to-noise ratio in the first pass. 
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New measurements of pulse-broadening times and predictions from the electron density models. 



PSR 


Ref a 


Period 


DM 


I 


b 


v 


PBFj 




PBF 2 




d , t c93 


a , tl & 2 00 . 






(ms) 


(pc cm~3) 


(°) 


C°) 


(MHz) 


T d (ms) 


St 


'j(-) 


St 


(ms) 


(ms) 


(!) 


(2) 


(31 


(4) 


(5) 


(6) 


(7) 


(S) 


(9) 


(10) 


(11) 


(12) 


(13) 


J1S4S + 0826 


1,400 


328.64 


90.8 


40.1 


4.6 


1175 


<0.2 




<0.2 




0.003 


0.001 


J1849+0127 


2,1400 


542. 11 


214.4 


33.9 


1.2 


430 


78±21 


0.9 


.*> 




13.4 


18.42 














1175 


6.5±2.1 


1.2 


3.3±1 


0.1 


0.161 


0.221 














1475 


<3.4 




<1.7 




0.059 


0.081 


J1850+0026 


2,1400 


1081.76 


197.2 


33.2 


0.5 


430 


9.6±2.4 


2.5 


4.8±1.1 


0.2 


10.56 


18.97 














1175 


<4.2 




<1.1 




0. 127 


0.228 


J1851 + 0118 


3,1400 


907.05 


413.0 


34. 1 


0.7 


1175 


<5 




<3 




1.027 


2.917 


J1851 + 0418 


1,300 


284.71 


112.0 


36.7 


2.1 


1175 


<3.5 




<1 




0.010 


0.008 


J18524-0031 


1,1400 


2180.06 


680.0 


33.5 


0.1 


1175 


495 ±25 


2.5 


271 ±5.7 


7.9 


4.846 


411.6 














1475 


225±14 


1.6 


127 ±3.7 


1.8 


1.782 


151.3 


J 1852 + 0305 


2,1400 


1326. 06 


3 .6 


3o.7 


1.3 


1175 


< 14 




. . .!> 




0.448 


0.728 


J 1853 + 0056 


2 1400 


275 56 


182 4 


33 9 


1 


117 r 


<2 




< 1 




089 


196 


J1853 + 0505 


3. 1400 


905.21 


273.6 


37.6 


2.0 


1175 


124 ±14 


0.5 


. . .*> 




0.243 


0. 193 














1475 


54 ±3 


1.4 


23±4 


1.8 


0.089 


0.071 


J1853+0545 


5,1400 


126.39 


198.7 


38.2 


2.3 


1175 


13.6 ±2 


0.4 


8.2 ±0.3 


2.5 


0.098 


0.049 














1475 


7.1 ±0.9 


0.3 


3 ±0.5 


0.8 


0.036 


0.018 














2380 


1.5 ±0.2 


0.1 


0.4 ±0.1 


1.7 


0.004 


0.002 


J1855 + 0422 


2,1400 


1677.98 


438.6 


37.2 


1.2 


1175 


27 ±3 


0.2 


16.6 ±2.5 


0.9 


1.004 


1.803 














1475 


< 12 




<4 




0.369 


0.663 


J1856+0113 


1,1400 


267.46 


96.7 


34.5 


-0.5 


1175 


<1 




<1 




0.006 


0.007 


J1856+0404 


2,1400 


420.22 


345.3 


37.1 


0.8 


1175 


9.5 ±4 


0.3 


4.8 ±1.6 


0.04 


0.709 


1.279 














1475 


6±3 


0.2 


2.8±1 


0.2 


0.261 


0.470 


J1857+0057 


1,300 


356.93 


83.0 


34.4 


-0.8 


430 


<5 




<1 




0.271 


0.229 


J1857+0210 


2,1400 


630.94 


783.2 


35.5 


-0.3 


1175 


13.4 ±3.6 


0.06 


6.1 ±0.65 


0.04 


9.004 


19.96 














1475 


<5.5 




<4 




3.311 


7.340 


J1857+0212 


1,1400 


415.80 


504.0 


35.5 


-0.2 


1175 


3.8 ±0.9 


0.1 


1.2 ±0.3 


1.8 


2.749 


6.135 














1475 


2.2±0.3 


0.8 


0.5±0.4 


0.7 


1.011 


2.256 


J1857+0526 


5,1400 


349.92 


468.3 


38.4 


1.2 


1175 


14.5 ±1.7 


0.4 


6 ±1 


0.6 


0.975 


1.700 














1475 


6.2±1.3 


0.04 


3±1 


0.03 


0.359 


0.625 


J1857+0809 


3,1400 


502.96 


284.2 


40.8 


2.5 


1175 


<3.5 




<2 




0.195 


0.061 


J1857+0943 




5.36 


13.3 




3.2 


430 










0.001 


0.001 


J1858+0215 


2,1400 


745.77 


702.0 


35.7 


-0.4 


1175 


38 ±3 


3.3 


22.5 ±6.5 


0.7 


6.267 


13.77 














1475 


18.4±4.7 


0.4 


11.7±4.8 


0.05 


2.304 


5.062 


J1858+0241 


3,1400 


4693.60 


341.7 


36.1 


-0.2 


1175 


<22 




<16 




0.758 


1.702 


J1900+0227 


2.1400 


374.24 


201.1 


36.1 


-0.8 


1175 


<4 




<2 




0.116 


0.178 


J1901 + 0156 


1,400 


288.22 


102.1 


35.7 


-1.2 


430 


3.5±1 


0.2 


0.7±0.3 


0.2 


0.619 


0.164 




























J1901 + 0331 


1,400 


655.48 


401.2 


37.2 


-0.5 


430 


60 ±3 


0.8 


44 ±4.2 


1.5 


107.9 


76.16 














1175 


<3 




<1 




1.295 


0.914 


J1901 + 0355 


3,1400 


554.80 


546.2 


37.5 


-0.3 


1175 


<4 




. . .*> 




2.97 


6.459 


J1901 + 0413 


2,1400 


2662.88 


367.0 


37.8 


-0.2 


430 


<558 




. . .*> 




79.5 


161.7 














1175 


< 13 




<3 




0.954 


1.940 


J1901 + 0716 


1,1400 


644.02 


252.8 


40.5 


1.2 


430 


10.1 ±2.4 


0.07 


8.5±3.4 


0.04 


17.12 


23.05 














1175 


<2.7 




<1 




0.205 


0.277 


J1901+1306 


1,400 


1830.72 


75.0 


45.7 


3.9 


1175 


. . d 




. . d 




0.001 


«1 


J1902 + 0556 


1,400 


746.60 


179.7 


39.4 


0.4 


430 


12 ±1.1 


1.1 


5.7 ±1.4 


0.3 


4.674 


10.05 




























J1902 + 0723 


1,400 


487.82 


105.0 


40.7 


1.0 


430 


<8 




<3 




0.455 


0.473 


J1903+0135 


1,400 


729.33 


246.4 


35.7 


-1.8 


430 


11.4 ±0.9 


3.1 


5.2 ±0.6 


3.5 


16.83 


13.73 














1175 


<3 




<1 




0.202 


0.165 


J 1903 + 060 1 


3,1400 


374.11 


398.0 


39.7 


0.2 


147 5 


2.7 ±0.5 


0.03 


1 ±0.4 


0.01 


1 405 


2.093 
















< 1 .7 




< 1 






0.77 


J1904+0004 


1,400 


139.53 


233.7 


34.4 


-2.8 


430 


3.1 ±1 


0.7 


4 ±1 


0.8 


10.78 


4.66 














1175 


<1.4 




<0.7 




0.129 


0.56 


J1904+0412 


2,1400 


71.09 


185.9 


38.1 


-0.9 


430 


■ ■ - d 




■ ■ d 




5.963 


10.24 


J1904+0800 


5,1400 


263.37 


438.3 


41.5 


0.9 


430 


. . . e 




. . . e 




61.02 


118.7 














1175 


3 ±0.3 


0.1 


1.2 ±0.12 


0.4 


0.723 


1.390 














1475 


<2 




< 1 




0.266 


0.511 


J1904+1011 


1,400 


1856.64 


135.0 


43.4 


1.9 


430 


. . b 




<4.4 




1.015 


0.937 


J1905+0400 


3,1400 


3.78 


25.8 


38.0 


-1.2 


430 


. . .d 




. . A 




0.005 


0.001 


J1905 + 0616 


2,1400 


989.64 


262.7 


40.1 


-0.2 


1175 


13.5 ±3.1 


0.04 


7.8 ±2.1 


0.2 


0.244 


0.514 














1475 


<5.8 




<1.5 




0.09 


0.189 


J1905 + 0709 


1, 1400 


648.05 


269.0 


40.8 


0.3 


430 


b 




41±10 


0.1 


20.95 


43.77 














1175 


7 ±4 


0.03 


3.2±1.6 


0.02 


0.251 


0.525 


J1906+0641 


1,1400 


267.29 


473.0 


40.5 


-0.2 


1175 


4.4 ±1.1 


0.2 


2.4±0.4 


0.7 


1.347 


111.1 














1475 


2.6±0.7 


0.05 


1.1±0.4 


0.12 


0.495 


40.85 


J1906+0912 


2,1400 


775.29 


260.5 


42.8 


1.0 


1175 


<8 




<4 




0.189 


0.278 


J1907+0534 


2,1400 


1138.31 


526.7 


39.7 


-0.9 


1175 


< 14 




<7 




1.351 


2.674 


J1907+0740 


2,1400 


574.65 


329.3 


41.5 


0.1 


430 


10.1±3.8 


0.03 


6.6 ±1.6 


0.03 


41.23 


85.56 














1 175 


< 2.3 




< 1 




0.495 


1.026 


J1907+0918 


1, 1400 


226. 12 


358.0 


43.0 


0.8 


1175 


<2.4 




<1.2 




0.471 


0. 766 


J1907+ 1247 


1,400 


827. 1 1 


257.0 






430 


b 




< 1 




7.631 


3.25 


J 1908 + 0457 








39 2 


l'*4 




. . .*> 






02 


















1175 


<3.6 




<1 = 8 




0.506 


0.69 


J1908+0500 


1,400 


291.03 


201.4 


39.3 


-1.4 


430 


4 ±0.7 


0.74 


2.1 ±0.25 


0.15 


7.703 


9.7 














1175 


<0.5 




<0.3 




0.092 


0.12 


J1908 + 0734 


1,400 


212.34 


11.1 


41.6 


-0.2 


1175 










«1 


«1 


J1908 + 0839 


2,1400 


185.38 


516.6 


42.5 


0.3 


1175 


5.6 ±1.3 


0.4 


3.5 ±0.8 


0.2 


1.131 


2.595 














1475 


2.6±0.6 


0.55 


1.4±0.2 


0.3 


0.416 


0.954 


J1908+0909 


2,1400 


336.53 


464.5 


43.0 


0.5 


1175 


4.9 ±0.9 


1.14 


2.4 ±0.5 


1 


0.757 


1.655 














1475 


<3 




. . b 




0.279 


0.609 


J1908 + 0916 


1,400 


830.31 


250.0 


43.1 


0.6 


430 


b 




12.7±2.1 


0.04 


12.87 


25.13 


J1909 + 0007 


1,400 


1016.96 


112.9 


35.0 


-3.9 


430 


<2 




<1 




1.006 


0.195 


J1909 + 0254 


1,400 


989.85 


172.1 


37.5 


-2.6 


430 


<2.7 




<1.1 




5.03 


1.995 


J1909 + 0616 


2,1400 


251.98 


348.6 


40.5 


-1.0 


430 


. . d 




. . d 




46.15 


69.79 














1175 


<5 




<2 




0.554 


0.837 


J1909 + 0912 


2,1400 


222.93 


421.4 


43.1 


0.3 


1175 


5.1±1.2 


0.04 


2.2 ±0.5 


0.03 


0.675 


1.681 














1475 


<2.2 




<1.1 




0.248 


0.618 


J1909+1102 


1,400 


283.65 


148.4 


44.7 


1.2 


430 


1.5 ±0.1 


2.4 


0.33 ±0.03 


5.2 


1.328 


1.076 



Bhat et al. 



Table 2 — Continued 



PSR Ref a Period DM I b u PBF, PBF 2 ^d,tc93 ^d,ne2001 

(ms) (pc cm- 3 ) (°) (°) (MHz) T d (m B ) / r T d (m B ) f r (ms) (ms) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) (10) (11) (12) (13) 



J1909 + 1450 


1,400 


996.04 


119.5 


48 


1 


2.9 


J1910+0358 


1,300 


2330.30 


78.8 


38 


6 


-2.3 


J 1 9 1 -j- 1 1 5 3 4 


2,1400 


452.83 


484.0 


40 





-1.6 


J1910 + 0714 


1,400 


2712.51 


124.1 


41 


5 


-0.8 


J1910 + 1231 


1,400 


1441.81 


274.4 


46 


2 


1.6 


.T1912 + 1036 


1,400 


409.38 


147.0 


44 


7 


0.3 


J1913+0446 


5,1400 


1616.09 


109.1 


39 


7 


-2.6 


,J1913-f-0832 m 


2, 1400 


134.40 


359.5 


43 





-0.9 


,T1913+0832 i 














J1913+0936 


1,400 


1242.06 


157.0 


43 


9 


-0.4 


noi q j_ i nnn 
.i iaio-|-i uuu 














J1913+1011 


2,1400 


35.91 


178.9 


44 


4 


-0.1 


J 19 13 + 1 145 


2 1400 


306 05 


630 7 


45 


g 


6 


J1913+1400 


1,400 


521.50 


144.4 


47 


8 


1.7 


J1914+1122 


1,400 


600.96 


80.0 


45 


6 


0.2 
















J1915+0738 


1,400 


1542.73 


39.0 


42 


4 


-1.7 


J1915+1009 


1,400 


404.56 


246.1 


44 


6 


-0.6 


.T1915 + 1606 


1,400 


59.06 


168.8 


49 


9 


2.2 


J1916+0844 


3,1400 


440.03 


339.0 


43 


6 


-1.1 


J1916 + 0951 


1,400 


270.26 


61.4 


44 


5 


-0.9 


J1916 + 1030 


1,400 


628.92 


387.0 


45 


1 


-0.6 


J1916 + 1312 


1,400 


281.86 


236.9 


47 


5 


0.7 


J1917+1353 


1,300 


194.62 


94.5 


48 


2 


0.8 


J191S | 1444 














J1918 + 1541 


1,400 


370.86 


13.0 


•19 


9 


1.4 


J1920+1110 


2,1400 


509.85 


181.1 


46 


1 


-1.2 


J1921 + 1419 


1,400 


618.14 


91.9 


49 





0.1 


.T1921 + 2003 


4,400 


760.70 


101.0 


54 


1 


2.8 


J1923+1706 


4,400 


547.23 


142.5 


51 


7 


1.0 


J1926 + 1434 


1,400 


1324.99 


205.0 


49 


8 


-0.8 


J1926 + 1928 


4,400 


1346.05 


445.0 


54 


1 


1.5 


J1927+1852 


4,400 


482.79 


254.0 


53 


7 


1.0 


J1927+1856 


4,400 


298.34 


90.0 


53 


8 


1.0 


J1929+1844 


4,400 


1220.38 


109.0 


53 


8 


0.5 


J1930+1316 


1,400 


760.05 


207.3 


49 


1 


-2.3 


J1931 + 1536 


4,400 


314.37 


140.0 


51 


3 


-1.4 


J1933 + 1304 


4,400 


928.38 


177.9 


49 


3 


-3.1 


J1935 + 1745 


4,400 


654.44 


214.6 


53 


6 


-1.2 


J1942+1743 


4,400 


696.30 


190.0 


54 


4 


-2.7 


J1944+1755 


4,400 


1996.78 


175.0 


54 


8 


-3.0 


J1945 + 1834 


4,400 


1068.80 


215.0 


55 


5 


-2.9 


J2027+2146 


4,400 


398.20 


96.8 


63 


5 


-9.5 



1175 <1 ■ ■ -° 0.016 0.013 



1175 


A 




A 




0.003 


0.012 


430 


< 4.7 




<5.7 




0. 14 


0.114 


1175 


12.5 ±3.8 


0.7 


5.4±2 


0.1 


0.676 


0.773 


1475 


<6.2 




<2.2 




0.249 


0.284 


430 


9.3±2.6 


0.01 


<1.1 




0.887 


1.628 


430 


24.6 ±6.4 


0.06 


14.7±1.8 


0.2 


10.59 


8.278 


1175 


<3 




<1.5 




0.127 


0.099 


430 


4.2±2.7 


0.04 


2.9±1.5 


0.03 


1.317 


1.933 


1175 


<3 




<2 




0.013 


0.023 


1175 


<2 




<1.4 




0.008 


0.005 


1175 


7.7 ±1.7 


0.7 


5.1 ±1.4 


0.1 


0.466 


0.734 




2.3±1 




. . A 








1175 


6.6 ±1.8 


0.12 


3 ±1.4 


0.2 


0.509 


0.905 


1475 


<2 




<1 




0.187 


0.333 


430 


<1.3 




<0.4 




1.874 


35.38 


1175 


1 1 . 1±4 


0.03 


5±3 


0.03 


0.556 


16.94 




<6 




<3 








430 


2±1 


0.3 


0.9±0.7 


0.4 


2.986 


107.9 


1175 


<0.4 




<0.2 




0.036 


1.294 


1175 


9.2 ±1 


0.3 


5.6±1.3 


0.1 


1.568 


1.204 


1475 


4.3±0.5 


0.2 


2.2±0.9 


0.01 


0.577 


0.443 


430 


3.3 ±0.6 


1.7 


<0.53 




0.706 


2. 126 


1175 


<2 




< 1 




0.009 


0.026 


1175 










0.001 


0.003 


1175 


<6 




<2 




0.396 


0.46 


430 


. . .c 




. . .c 




0.012 


0.002 


430 


15.4 ±1.3 


6.9 


11.1 ±1 


0.9 


10.20 


20.27 


1175 


< 1 ■ 1 




<0.4 




0.122 


0.243 


1175 


0.33±0.1 


0.4 


<0.07 




0.007 


0.021 
















1475 


3.6 ±0.9 


0.1 


1.1 ±0.4 


0.3 


0.131 


0.158 


430 


1±0.4 


0.6 


<0.6 




0.037 


0.025 


1175 


9.2±2.9 


0.02 


4.8 ±1.5 


0.02 


0.4 


0.714 


1475 


<4 




<2 




0.147 


0.262 


430 


2.8 ±0.1 


5.7 


0.9 ±0.4 


2.4 


5.958 


12.93 


1175 


<1 




<0.4 




0.071 


0.155 


1175 


1.2±0.4 


0.3 


0.4±0.1 


0.4 


0.001 


0.011 


430 


. c 




. c 




0.006 


0.001 


1175 


. c 




. c 




«1 


<£1 


1175 


14 ±3.3 


0.02 


7.4±1.4 


0.01 


0.032 


0.06 


1475 


6.3±2.1 


0.4 


2.9±1.7 


0.5 


0.012 


0.022 


1175 


6±4 


0.4 


3.2±1.9 


0.9 


0.002 


0.01 


1475 


4±3 


0.3 


2.9±2.2 


0.05 


0.001 


0.004 


430 


. . .!> 




<1.2 




0.096 


0.07 


430 


. . A 




<1 




0.406 


1.23 


430 


8.4±1.3 


4.2 


<7.4 




1.934 


5.43 


430 


. . A 




34.5 ±6.8 


0.4 


11.46 


3.46 


430 


. . b 




5.8 ±1.3 


0.5 


2.754 


1.97 


1175 


. . A 




<1 




0.033 


0.024 


430 


. . A 




<0.2 




0.09 


0.09 


1175 


. . b 




<1.1 




0.002 


0.003 


430 


. . A 




4.1±2.4 


0.1 


1.526 


1.89 


1175 


. . A 




<1.2 




0.018 


0.023 


430 


. . A 




0.9 ±0.2 


0.01 


0.33 


0.87 


430 


A 




<0.25 




0.660 


0.334 


430 


15.9 ±1.5 


1.2 


9.2±3.3 


0.02 


0.751 


2.72 


1175 


<1.1 




<1.2 




0.008 


0.033 


430 






5±1.5 


0.5 


0.403 


0.29 


430 


. . b 




<5 




0.308 


0.215 


430 


. . A 




3.3 ±1.3 


0.2 


0.7 


0.524 


430 


<0.4 




<0.4 




0.091 


0.041 



a Reference(s) to pulsar parameters, followed by the frequency band (in MHz) of the survey that discovered the pulsar. (1) ATNF pulsar catalog, 
available at http://www.atnf.csiro.au/rescarch/pulsar/psrcat; (2) Morris et al. (2002); (3) Hobbs ct al. (2004), in preparation; (4) Lorimcr ct al. 
(2002): (5) Kramer ct al. (2003). Note that pulsars in (2), (3), and (5) (Parkes multibeam survey discoveries), as well as in (4), are also available in 

Ci). 

°The PBF yields unphysical residuals from the deconvolution for any realistic value of t ^ (sec text). 
c Signal-to-noisc ratio is too small to allow a meaningful fit to the PBF. 
d T d is negligibly small. 

e The pulsar is not detected at this frequency. 

m,1 Main pulse (m) and inter pulse (i) of the pulsar (see also Fig. 1). 



Multifrequency Observations of Pulse Broadening 



Table 3 

Frequency scaling indices from measurements of pulse-broadening times. 



PBFi deconvolution PBF2 deconvolution 



PSR 


(MHz) 


V2 

(MHz) 






a 1 






(h 


&2 


02 


(1) 


(2) 


(3) 






(4) 






(5) 


(6) 


(7) 


J1849+0127 


430 


1175 


2 


.5 


± 


0. 


1 


a 






J1852+0031 


1175 


1475 


3 


.5 


± 


0. 


1 


4.7 ± 0.4 


3.3 ± 0.1 


5 ± 0.2 


J1853+0505 


1175 


1475 


3 


.7 


± 


0. 


2 


4.4 ± 0.5 






J1853+0545 


1175 


1475 


2 


.8 


± 


0. 


3 


6.8 ± 1.5 


4.4 ± 0.3 


3.7 ± 0.5 




1475 


2380 


3 


.2 


± 


0. 


1 


5.3 ± 0.5 


4.2 ± 0.2 


3.8 ± 0.3 




1175 


2380 


3 


.1 


± 


0. 


1 


5.7 ± 0.3 


4.3 ± 0.1 


3.8 ± 0.1 


J1856+0404 


1175 


1475 


2 


.0 


± 


0. 


.8 


a 


2.4 ± 0.7 




J1857+0212 


1175 


1475 


2 


.4 


± 


0. 


4 


a 


3.9 ± 1.4 


4.1 ± 3.2 


J1857+0526 


1175 


1475 


3 


.7 


± 


0. 


4 


4.3 ± 1 


3 ± 0.6 


6.1 ± 2.6 


J1858+0215 


1175 


1475 


3 


.2 


± 


0. 


.4 


5.4 ± 1.6 


2.9 ± 0.8 


6.6 ± 3.8 


J1906+0641 


1175 


1475 


2 


.3 


± 


0. 


5 


a 


3.4 ± 0.6 


4.8 ± 2 


J1908+0839 


1175 


1475 


3 


.4 


± 


0. 


5 


4.9 ± 1.7 


4.0 ± 0.4 


4 ± 1 


J1913+0832 


1175 


1475 


5 


.3 


± 


0. 


.8 


3.2 ± 1.1 






J1913+1145 


1175 


1475 


3 


.4 


± 


0. 


3 


5 ± 0.8 


4.1 ± 0.8 


3.9 ± 1.6 


J1916+0844 


1175 


1475 


3 


.3 


± 


0. 


A 


5 ± 1.5 


5.7 ± 0.6 


3.1 ± 0.7 


J1920+1110 


1175 


1475 


3 


.5 


± 


0. 


1 


4.7 ± 1.9 


4.1 ± 1 


3.9 ± 2.1 


J1921+1419 


1175 


1475 


1 


.8 


± 


1. 


2 


a 


0.4 ± 0.4 





a Implicd values for (3 are unphysically large. 
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Table 4 

Estimates of scattering measures and constraints on the properties of clumps. 



PSR 


DM 


D 


I 


b 


Freq. 


SM meas 


SM„ e2 001 


(5SM 


F C (5DM) 2 




(pc cm °) 


(kpc) 


(deg) 


(deg) 


(MHz) 


(kpc m M ' 6 ) 


(kpc m M ' 6 ) 


(kpc m /U/J ) 


(log) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


J1849+0127 


214.4 


5.49 


33.95 


1.20 


430 


0.21 


0.11 


0.093 


-4.04 


J1853+0546 


197.2 


5.37 


38.24 


2.28 


1175 


3.58 


0.033 


3.54 


-2.45 












1475 


4.79 




4.76 


-2.32 












2380 


7.67 




7.64 


-2.11 


J1855+0422 


438.6 


8.39 


37.22 


1.20 


1175 


4.36 


0.46 


3.90 


-2.60 


J1856+0404 


345.3 


6.85 


37.07 


0.84 


1175 


2.15 


0.41 


1.74 


-2.86 


J1857+0526 


468.3 


8.92 


38.40 


1.24 


1175 


2.47 


0.41 


2.05 


-2.90 












1475 


2.33 




1.92 


-2.93 


J1858+0215 


702.0 


10.02 


35.68 


-0.43 


1175 


4.98 


2.14 


2.84 


-2.81 


J1901+0413 


367.0 


6.81 


37.77 


-0.20 


430 


4.32 


0.58 


3.75 


-2.53 


J1903+0609 


398.0 


7.21 


39.72 


0.24 


1175 


0.73 


0.59 


0.14 


-3.98 


J1904+0802 


438.3 


8.67 


41.51 


0.89 


1175 


0.68 


0.36 


0.32 


-3.69 


ti 9054-O61 r 


969 7 


5 76 

O . 1 U 


40 O^i 

L ±\) .\J'J 


-0 1 5 


1 1 75 


3 35 


0.22 


3 1 3 


-9 53 


J1908+0839 


516.6 


9.35 


42.51 


0.29 


1175 


1.07 


0.57 


0.51 


-3.53 


J1908+0909 


464.5 


8.96 


42.96 


0.52 


1175 


0.99 


0.40 


0.60 


-3.44 


J1909+0912 


421.4 


8.27 


43.11 


0.32 


1175 


1.11 


0.44 


0.67 


-3.36 


J1910+0534 


484.0 


10.40 


40.00 


-1.57 


1175 


1.92 


0.19 


1.73 


-3.04 


T1 Ql 94-0828 


359 5 

OOiJ . o 


7 74 


42.81 


-0 67 


1 1 75 


80 


98 


53 


-3 43 


J1913+0832 


359.5 


7.93 


42.98 


-0.86 


1175 


1.60 


0.23 


1.38 


-3.03 


J1913+1145 


630.7 


14.56 


45.83 


0.63 


1175 


1.12 


0.21 


0.92 


-3.47 


J1920+1110 


181.1 


5.60 


46.11 


-1.16 


1175 


3.53 


0.038 


3.49 


-2.47 


J1852+0031 


680.0 


7.91 


33.46 


0.11 


1175 


51.65 


44.29 


7.36 


-2.30 












1475 


61.52 




17.23 


-1.93 


T1 009-1-0^^ 


1 70 7 


4 70 


30 41 


"^fi 
u.ou 


zno 


OQ 


078 
U.u ( o 


01 3 


-4 84 


J1909+1102 


148.4 


4.17 


44.74 


1.17 


430 


0.017 


0.013 


0.0038 


-5.31 


J1910+0714 


124.1 


4.05 


41.48 


-0.80 


430 


0.082 


0.019 


0.063 


-4.07 


J1910+1231 


274.4 


7.73 


46.17 


1.63 


430 


0.11 


0.044 


0.065 


-4.34 


J1913+1400 


144.4 


5.12 


47.82 


1.67 


430 


0.029 


0.020 


0.0089 


-5.03 


J1915+1606 


168.8 


5.91 


49.91 


2.22 


1175 


0.15 


0.015 


0.13 


-3.91 


J1916+1030 


387.0 


8.59 


45.06 


-0.60 


1175 


1.74 


0.21 


1.54 


-3.01 


J1926+1434 


205.0 


6.44 


49.80 


-0.84 


430 


0.052 


0.036 


0.016 


-4.88 


J1853+0505 


273.6 


6.49 


37.64 


1.97 


1175 


19.20 


0.088 


19.11 


-1.80 












1475 


22.17 




22.08 


-1.73 


J1915+0856 


339.0 


7.96 


43.56 


-1.11 


1175 


1.60 


0.14 


1.46 


-3.00 












1475 


1.93 




1.78 


-2.92 


J1935+1745 


214.6 


6.93 


53.63 


-1.21 


430 


0.082 


0.019 


0.06 


-4.30 
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=1849+0127 
== 542 ms 





430 MH: 
DM=214.4 



pms+sm 

== 542 ms 



1175 MHz 
DM=214.4 



*1S49+012? 
== 542 ms 


' 1475 MH: 
DM=214.4 







PiM+bW 23*6 mhz| 

== 542 ms DM=214.4 




=1850+0026 


430 MHi 


==1082 ms 


DM= 197.2 









=!SSd+bb2B 

==1082 ms 

' 


' 1175 MHi 
DM= 197.2 





=1850+0026 


, ' 1475 MH: 


==1082 ms 


DM=197.2 









=1850+0026 , ' 2380 MHz 
==1082 ms DM= 197.2 



prnmwr '1^75™ PTgsT+orre — r -~Tmrm 

== 907 ms I DM=413.C == 907 ms DM=413.C 




J185T+0418 
== 285 



430 MH: 
DM=112.C 





J1851+0418 , ' 1475 MH; 
== 285Ans f[ DM=112.C 





J 1852+0031 
==2180 ms 


„ ' 1175 MH: 
Ht DM=680.C 







J1852+0031 
==2180 ms 


' 1475 MHi 
A DM=680.C 

X 





J1 852+0031 
==2180 ms 


, ' 2380 MH: 
f| DM=680.G 







==1326 ms 





1175 MHi 
DM=315.6 




==1326 ms f, DM=315.E ==1326 ms 







=1853+0056 „ ' 1475 MHi 
== 276 ms I DM= 182.4 




=1853+0505 , 


1475 MH: 


== 905 ms 


DM=273.e 





PTI53+5B56" , 

== 276 ms 




2380 MHi 
DM= 182.4 




= 1853+0505 I 2380 MHi 
== 905 ms J| DM=273.e 





=^3+0545 


i ' 1175 MHi 


== 126 ms 


A DM=197.2 







il'dS3+054B 


' 1475 MH: 


'= 126 ms 


DM=197.2 







== 126 ms 


' 2380 MHi 
i DM= 197.2 





=1855+0422 
==1678 ms 




=1855+0422 
==1678 ms 


' 1475 MHi 
DM=438.C 







=1855+0422 
==1678 ms 


' 2380 MH: 
DM=438.e 
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J1856+0113 
== 267 ms 


' 1175 MHz 
DM= 96.7 







J1 857+0057 
== 357 ms 


, 430 MHi 
i| DM= 83.0 







=1856+0404 . ' 1175 MHi 
3 = 420 ms , DM=345.: 




J 1857+0057 
== 357 ms / 


' 1 1 75 MHi 
\ DM= 83.C 







J 1857+0057 
3 = 357 ms I 


' im m 

\ DM= 83.C 








J1 857+621 2 
== 416 ms 


' 1475 MHi 
DM=504.C 






J1857+0212 


' 2380 MHi 


== 416 ms 


DM=504.C 







^&7+052B 




== 350 ms 









=1857+0809 



430 MHz 
284.2 




=1858+0241 
=4694 ms 



430 MHi 
DM=341.7 




PTisT+BloTT - ' n75~'MhTa 

== 503 ms f, DM=284.2 




=18iS7+b526 


uww 


== 350 ms 


DM=468.: 









=1858+0241 
==4694 ms 


' 1175 MHi 
DM=341.7 







=1857+0809 ' 1475 MHi 
== 503 ms I DM=284.2 




== 746 ms 1 DM=702.C == 746 ms I DM=702.C 





== 374 ms 

.1 i.lu Lil.JlLlnriJ IK JjJ 


' 1175 MHi 
DM=201.1 





=1900+0227 I ' 1475 MHi 
== 374 ms DM=201.1 




=1857+0526 
== 350 ms 

: I 


' 2380 MHi 
DM=468.3 








Multifrequency Observations of Pulse Broadening 



J1 901+01 56 
3 = 288 ms 

„ j 


430 MHz 
DM= 102.1 





J1 901+0331 
3 = 655 ms 

J 


. ' 430 MHz 
\ DM=401.2 






J1 902+0556 


430 MHz 


== 747 ms 


DM= 179.7 









J1 902+0723 , ' 430 MH: 
=>= 488 ms DM=105.C 




J1 903+01 35 
== 729 ms 

_.. . 


430 MHz 

H 





I190T+0156 . ' 1175MHz 
s = 288 ms DM= 102.1 




J190T+0331 
3 = 655 ms 

I 


' 1175MHz 
DM-401.2 

\ 





PT33T+6555" 1 , if75 , MH- 

>- 555 ms I DM=546.2| 



3 =2663 ms 


' 1175MHz 
^ DM=367.0 







11901+0716 
3 = 644 ms 


' 1175MHz 
DM=252.£ 





J 1902+0556 
»= 747 ms 


' 1175MHz 
DM= 179.7 






J1903+0135 
3 = 729 ms 



1175 MHi 
DM=246.4 



="1 903+0601 . ' 1175 MHz 
== 374 ms DM=398.C 



J1901+0156 
?- 288 ms 



1475 MHz 
DM=102.1 




J1901+0331 
== 655 ms 

' _ i 


■ 1475 MHz 
DM=401 .2 





=W<M+635S 


14/5 MHj 


3 = 555 ms 


DM=546.2 







3T35T+54TT 

=>=2663 ms 



1475 MHz 
DM=367.C 



J1901+0716 

n 


■ 1475 MHz 
DM=252.E 





J 1902+0556 
°= 747 ms 

„. .. 


1175 MHz 
DM=179.7 





11902+0723 
488 ms 



1475 MHi 
DM=105.C 



J1903+0135 
3 = 729 ms 




1475 MHz 
DM=246.4 

I 





?1 903+0601 ," 1475 MHj 
== 374 ms , DM=398.C 





J1 901+01 56 
== 288 ms 



2380 MHz 
DM= 102.1 





^dl+0413 


' 2380 MHz 


==2663 ms 


DM=367.C 




L-; 



J1 901+071 6 


' 2380 MHz 


3 = 644 ms 


DM=252.£ 







J1903+0135 
== 729 ms 


2380 MHz 
DM=246.4 

I 





=1903+0601 


' 2380 MHz 


== 374 ms 


DM=398.C 







J1 904+0004 . 430 MHz 
== 140 ms h DM=233.7 



J 1904+0004 ' 1175 MHz 
3 = 140 ms DM=233.7 



J 1904+0004 ' 1475 MHz 
== 140 ms A DM=233.7 



J1 904+0004 [ ' 2380 MHz 
°= 140 ms , DM=233.7 
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J1 908+0457 
== 847 ms 



430 MHz 
DM=360.C 




J1 908+0500 
== 291 ms 



430 MHz 
DM=201.4 



J1 908+0916 
== 830 ms 



430 MHz) 
DM=250.C 




J1 909+0007 
==1017 ms 


430 MHz 
DM= 112.9 

L. - 





J1 909+0254 
== 990 ms 


430 MHz 
DM= 172.1 

\ 





=1909+0616 , 430 MH: 
== 252 ms i A DM=348.e 




11908+0457 


' 1175MH2 


== 847 ms i 


DM=360.C 









J 1908+0500 


■ 1175MHz 
DM=201.4 





J 1908+0734 , ' 1175 MH; 
== 212 ms 




=1908+0839 . ' 1175 MHz 
== 185 ms \ DM=516.6 





=1908+0909 
== 337 ms 



1175 MHz 
DM=464.S 



11908+0916 . ' 1175 MHzl 
== 830 ms DM=250.C 




J 1909+0007 ' 1175 MHz 
==1017 ms DM=112.9 



J 1909+0254 
== 990 ms 



1175 MHz 
DM= 172.1 




J 1908+0457 


' 1475 MHz 


== 847 ms 


DM=360.C 









J 1908+0500 
== 291 ms 


' 1475 MHz 
DM=201 .4 





J 1908+0734 


1 475 MHz 


== 212 ms 


DM= 11.1 


Z) 







=>1 9(08+0839 
== 185 ms 


' 1475 MHz 
DM=516.e 







=1908+0909 


' 1475 MHz 


== 337 ms 


DM=464.J 







J 1908+09 16 . ' 1475 MHzl 
== 830 ms I DM=250.C 




11909+0007 ' 1475 MHz 
1017 ms DM=112.9 



J 1909+0254 
== 990 ms 


1475 MHz 
DM=172.1 





=1909+0616 ,, ' 1475 MHzl 
== 252 msfi DM=348.e 




=1909+0912 L 
== 223 ms f 


' 1475 MHz 
DM=421 A 





J1 908+0457 


' 2380 MHz 


== 847 ms 


DM=360.C 







J1 908+0500 
== 291 ms 


' 2380 MHz 
DM=201.4 






J 1909+0254 , 2380 MHz 
== 990 ms DM= 172.1 




*«d9+b9i2 

== 223 ms 


' 2380 MHz 
DM=421.4 




mm 
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J1913+1400 
== 521 ms 


430 MHz 
DM= 144.4 

\ 





J1913+1400 
== 521 ms 

"_ \ 


' 1175 MHz 
DM= 144.4 

V . _ . .... .. .. 





J1913+1400 
== 521 ms 

J 


' 1475 MHz 
DM=144.4 





J1913+1400 I ' 2380 MHz 



== 521 ms 



DM= 144.4 



11914+1122 
3 = 601 ms 


, ' 1 1 75 MHz 
DM= 80.C 







J1914+1122 
== 601 ms 


■ 1475 MHz 
DM= 80.C 

| 







J1914+1122 ' 2380 MHz 
=>= 601 ms DM= 80.C 




J1915+0738 
==1543 ms 


430 MHi 
DM= 39.0 


IliiiLnLjlHinij 

ifipiwBPOIj 





J1 91 5+0738 
'=1543 ms 


' 1175MHz 
DM= 39.C 





J1 91 5+0738 
==1543 ms 


1 475 MHz 
DM= 39.C 





PT3T5+61S8" ",' " 1175MHz| 
== 343 ms ft DM=369.1 




P1SH5+6839 ' "1475 MHil 
== 343 ms DM=369.1 





== 343 ms J DM=369.1 






=W<Si6+b844 
== 440 ms 


' 1175MHz 
DM=339.C 





=1916+0844 
3 = 440 ms 


' 1475 MHz 
DM=339.C 

I 






J1915+1009 


430 MHz 


== 405 ms 


\ DM=246.1 







J1915+1009 
== 405 ms 


if 75 MHz 
DM=246.1 

.._ 





J1915+1009 
°= 405 ms 

... . 


T475 MHz 
DM=246.1 





J1915+1009 
== 405 ms 


' 2380 MHz 
DM=246.1 





11915+1606 . ' 1175 MHi 
>= 59 ms i DM=168.£ 





11916+0951 
== 270 ms 


430 MHz 
DM= 61.4 





11916+0951 
== 270 ms 

... _ . i 


, 1175 MHz 
DM= 61.4 





J1916+1030 " 1175 MHz 
'= 629 ms 1 DM=387.C 




J1916+0951 
== 270 ms 


1475 MHz 
DM= 61.4 






J1916+0951 
== 270 ms 


2380 MHz 
DM= 61.4 







J1916+1312 

'"I 


430 MHz 





J1916+1312 
== 282 ms 


' 1175 MHz 
DM=236.S 





J1916+1312 
== 282 ms 


' 1475 MHz 
DM=236.S 





J1916+1312 
== 282 ms 


' 2380 MHz 
DM=236.S 
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J1917+1353 
== 195 ms 

J 


' 1175MHz 
DM= 94.5 





J1918+1444 
==1181 ms 


430 MHz 
DM= 30.C 





J1926+1434 


' 430 MHj 


==1325 ms 


DM=205.0 








J'1 930+1316 
== 760 ms 



J1918+1444 
==1181 ms 




' 1175 MHz 
DM= 30.0 






J1926+1434 
==1325 ms 


' 1175 MHz 
DM=205.0 

A ; 





1175 MHz 
DM=207.J 



J1917+1353 
== 195 ms 


1475 MHz 

1"" 





J1918+1444 
==1181 ms 


' i 475 MHz 
DM= 30.C 





=1920+1110 
== 510 ms 


' 1475 MHz 

I ™ : 







J1921 + 1419 


f 475 MHz 


s = 618 ms 


DM= 91 .S 







J 1926+ 1434 
==1325 ms 


' 1475 MHz 
DM=205.C 




J1930+1316 
== 760 ms 


' 1475 MHz 
DM=207.2 



J1917+1353 
== 195 ms 


' 2380 MHz 
DM= 94.E 







J1918+1444 


' 2380 MHz 


==1181 ms 


DM= 30.C 








J1926+1434 
'=1325 ms 



2380 MHz 
DM=205.C 



Fig. 1. — Integrated pulse profiles of pulsars from Arecibo observations at 430, 1175, 1475 and 2380 MHz. Data at 430 MHz were taken with 
the PSPM, and those at higher frequencies were taken with the WAPP. All profiles are plotted with a pulse phase resolution of 2 milli-periods, 
where consistent with the data acquisition time resolution (Table 1). The highest point in the profile is placed at phase 0.5. The pulsar ID, 
period and the dispersion measure are indicated at the top of each panel, along with the center frequency of observation. Objects with labels 
(top left) starting with 'P' refer to new discoveries from the Parkcs multibeam survey, and those with 'J' are previously known pulsars. 
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P1913 + 0832(ip) 



P1913+0832(mpj , 


1 175 MHz 


1.5 




Li 


u 


1 

0.5 

















50 


100 




J 1 91 5 + 1 009 


430 MHz 


1.5 








1 

0.5 










1 1 75 MHz 




P1913+1 145 



J 1 9 1 6+ 1 030 - 


1 175 MHz 




u 


u 







1 1 75 MHz 



P1920+1 110 



100 200 



1 175 MHz 




Fig. 2. — Examples of the intrinsic pulse shapes (light solid [red] curves with highest peaks) and the best fit PBFs (solid [blue] curves rising 
from zero at left of each panel) obtained by application of the CLEAN method; the PBF is assumed to be a simple one-sided exponential 
(PBFi, appropriate to a thin slab scattering geometry). The amplitudes of both the PBFs and the measured profiles (heavy solid [green] 
curves) are normalized to unity, and the areas under the intrinsic and measured pulse profiles are identical. 
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J1852 + 0031 


■ 1 1 75 MHz 


" ■ •■• " "V' 





J1852+0031 


1475 MHz 




= 1 853 + 0545 


1 175 MHz 






5 














P1 853+0545 
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Fig. 3. — Similar to Figure 2, except that the pulse broadening function employed by CLEAN has a more rounded shape (PBF2, due to 
uniform scattering medium between the pulsar and the Earth). 
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10 100 1000 



DM (pc cm ) 

Fig. 4. — Measurements of pulse-broadening times plotted against dispersion measures. The new measurements are shown as filled circles. 
The open circles with crosses (DM < 200 pc cm -3 ) are derived from the measurements of decorrelation bandwidths, while the open circles 
are published measurements. The solid curve represents the best fit model for the empirical relation between T,j and DM, the frequency- 
independent coefficients for which are only slightly different from those obtained by Cordes & Lazio (2002b) based on the published data 
alone (see § 5.2 for details). 
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Fig. 5. — Analysis of the frequency dependence of pulse broadening that takes into account an inner scale for the wavenumber spectrum 
of electron-density irregularities. (Top panel) Plot of 8a, the difference in exponent in the relation r d oc v~ a above and below a break point 
defined by the composite quantity Td tCIOaa v 2 / D. We calculate the best fit values of a for data points above and below the break point and 



calculate 8a as a function of r d , 



2 /D. The units of cross f 2 /D are (ms GHz 2 kpc 1 ). (Bottom panel) \ 2 f° r the fit as a function of 



T d cross 17 /D, defined here as the sum of the squares of data— model (see text). 
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Fig. 6. — Measurements of pulse-broadening times plotted against the predictions from the new electron density model NE2001 (Cordes 
& Lazio 2002a). The filled circles are the published measurements. The new measurements from our observations are shown as open circles. 
All measurements arc scaled to a common frequency of 1 GHz using Td <x v~ 4,4 . The dashed line is of unity slope. As evident from the 
figure, a significant number of both the published and new measurements are well above the dashed line, which implies that the model tends 
to underestimate the degree of scattering toward many lines of sight. 
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Fig. 7. — Measurements of pulse-broadening times plotted against the predictions from the electron density model TC93 (Taylor & Cordes 
1993). The filled circles are the published measurements, and the new measurements from our observations are shown as open circles. All 
measurements are scaled to a common frequency of 1 GHz using Td c< v~ iA . The dashed line is of unity slope. As for Figure 6, the model 
tends to underestimate the degree of scattering toward many lines of sight. 
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Fig. 8. — Estimates of scattering measure (SM) derived from all pulse-broadening data available. The size of the symbol is proportional to 
log(SM). Pulsar positions are projected onto the Galactic plane; filled circles represent the published data, and the new measurements from 
our observations are shown as open circles. The spiral arm locations are adopted from the NE2001 model of Cordes & Lazio (2002a, b). 
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Fig. 9. — Similar to the plot in Figure 8, except that the quantity plotted is the departure of the measured pulse-broadening time (r^) 
from the prediction of the NE2001 model (t^ ne200l); the size of the symbol is proportional to the absolute value of log(T ( j/T t j ne 200l)- As for 
Figure 8, the filled circles represent the published data, while the open circles are the new measurements. 
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Fig. 10. — Estimates of .F C (<5DM) 2 for the clumps of enhanced scattering, derived from the excess scattering measures (Table 4), are plotted 
against the respective pulsar dispersion measures. The results are for a clump size of ~ 10 pc and a volume number density ~ 1 kpc - 3 for the 
clumps. For a fluctuation parameter of F c = 10, these results imply excess DM within the range 7 X 10 — 4 -4 X 10 — 2 pc cm -3 . 
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Fig. 11. — Measurements of frequency scaling index (cti) against the respective DMs. The results for low-DM objects (Cordes, Weisberg, 
& Boriakoff 1985) are derived from measurements of decorrelation bandwidths. For PSR J1852+0031, the only object common between 
our sample and that of Lohmer et al. (2001), estimates of a are consistent within measurement errors. The dashed line corresponds to the 
Kolmogorov scaling index. 



